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Dynamical chiral symmetry breaking in the DLCQ method 
is investigated in detail using a chiral Yukawa model closely 
related to the Nambu-Jona-Lasinio model. By classically 
solving three constraints characteristic of the light-front for- 
malism, we show that the chiral transformation defined on 
the light front is equivalent to the usual one when bare 
mass is absent. A quantum analysis demonstrates that a 
nonperturbative mean-field solution to the "zero-mode con- 
straint" for a scalar boson a can develop a nonzero conden- 
sate (a) = — ■j^($\E') 7^ while a perturbative solution cannot. 
This description is due to our identification of the "zero-mode 
constraint" with the gap equation. The mean-field calculation 
clarifies unusual chiral transformation properties of fermionic 
field, which resolves a seemingly inconsistency between trivi- 
ality of the null-plane chiral charge Qs F |0) = and nonzero 
condensate (4'^/) ^ 0. We also calculate masses of scalar and 
pseudoscalar bosons for both symmetric and broken phases, 
and eventually derive the PCAC relation and nonconservation 
of Q\ f in the broken phase. 

PACS number(s): 11.30.Rd, 11.30.Qc, 11.15.Pg, 12.40.-y 



I. INTRODUCTION 

Chiral symmetry breaking is undoubtedly one of 
the most important concepts for understanding hadron 
physics in low energy region jlj . The smallness of n and 
K masses is beautifully explained if one identifies them 
with the Nambu-Goldstone (NG) bosons associated with 
the chiral symmetry breaking. An important aspect of 
this phenomenon is dynamical formation of NG bosons 
as bound states of quarks and gluons in the strong cou- 
pling region. However, its complete demonstration in 
QCD is not reached yet because of the difficulties in 
describing bound states in a nonperturbative and rela- 
tivistic manner. Instead, many people have been inves- 
tigating much simpler effective models of QCD. Among 
them, the Nambu-Jona-Lasinio (NJL) model || is the 
most deeply and thoroughly understood. The NJL model 
is a 3+1 dimensional four- Fermi theory and reproduces 
various properties of hadrons concerning the chiral sym- 
metry breaking despite some undesirable features such as 
nonrenormalizability and lack of confinement [^) . Nowa- 
days the model plays a role of laboratory in which we 
can test new ideas proposed for nonperturbative study of 



low energy QCD. Therefore the NJL model is the most 
appropriate model in which we can check whether the 
light-front (LF) quantization can be applied to the dy- 
namical chiral symmetry breaking. The LF quantization 
is a newly revamped nonperturbative method for solving 
relativistic bound states in quantum field theory 0. 

Let us explain why the chiral symmetry breaking be- 
comes a special issue in the LF formalism. The rea- 
son is twofold: the first is apparent contradiction be- 
tween a nontrivial vacuum and a LF "trivial" vacuum, 
and the second is peculiarity of LF chiral transforma- 
tion. To resolve these problems is our primary purpose 
in the present paper. One of the remarkable merits of 
the LF quantization is that the Fock vacuum defined by 
a free theory is also the vacuum of the full theory. Many 
technical advantages such as exact Fock state expansion 
arise from this fact. On the contrary, the conventional 
formulation says that the chiral symmetry breaking is es- 
sentially a physics of finding another vacuum that breaks 
the chiral symmetry but is energetically favored. Such 
"vacuum physics" is thought to be very important for 
understanding nonperturbative phenomena in low energy 
region. Therefore to apply the LF formalism to QCD 
necessarily entails a problem how to realize such "vac- 
uum physics" within a framework with a trivial vacuum. 
For the purpose of understanding this problem, there 
are considerable efforts || to describe the spontaneous 
symmetry breaking in a simple scalar model (A</>f +1 ). 
They succeeded in obtaining the critical coupling which 
is consistent with the conventional results. The key is to 
solve a constraint equation for the longitudinal zero mode 
("zero-mode constraint") which appears in the DLCQ 
(Discretized Light-Cone Quantization) method 0. A 
nonzero condensate is realized as a nonperturbative so- 
lution of the zero- mode constraint. We will discuss this 
method in more detail later. 

Compared with such extensive studies, only little 
is known about the dynamical symmetry breaking in 
fermionic systems. Especially there have been only few 
attempts about the NJL model on the LF ]7|-[l0|] . At first 
glance, it seems not possible to follow the same route as in 
the scalar models because we do not have bosonic fields 
as fundamental degrees of freedom in the NJL model. 
However, we can apply the same idea to the dynamical 
symmetry breaking if one introduces bosonic auxiliary 
fields to the fermion bilinears and raises them to dynami- 
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cal variables by adding their kinetic terms. Of course the 
original fermionic model is reproduced as an infinitely 
heavy mass limit of the bosonic fields. According to this 
idea, we succeeded in describing the dynamical symmetry 
breaking (discrete chiral symmetry) in the 1+1 dimen- 
sional four- Fermi theory (the Gross-Neveu model) |pl| . 
The present paper is a generalization of this preliminary 
work which discussed only discrete chiral symmetry. We 
consider a kind of Yukawa model with continuous chiral 
symmetry, which is obtained from the NJL model using 
the above technique. We work within the DLCQ method 
so that we can formulate the problem from the viewpoint 
of the zero-mode constraints. It should be commented 
however that it is possible to discuss the dynamical sym- 
metry breaking even without introducing auxiliary fields. 
In Ref . H , one of the authors insisted the importance of 
a "fermionic constraint" which is again unique to the LF 
formulation and has very complicated structure due to 
the four-Fermi interaction. (Another merit of including 
scalar fields is a quite simplification of the fermionic con- 
straint.) More detailed analysis in this direction will be 
reported in the next paper 

One more point to be discussed is the unusual behavior 
of chiral transformation on the LF. In the LF formula- 
tion, a half degree of freedom of the fermion is a de- 
pendent variable to be represented by other independent 
variables. Therefore, chiral transformation should be im- 
posed only on the independent component of the fermion 
p3[ . It is not clear in interacting models whether the LF 
chiral transformation is equivalent to the usual one or 
not. 

The paper is organized as follows. In the next section, 
we define the chiral Yukawa model which is closely re- 
lated to the NJL model and introduce our framework, 
the DLCQ method. The classical aspects of the model 
is discussed in Sec. III. Here, we see that there are 
three constraints (i.e. two zero-mode constraints and one 
fermionic constraint). We also show peculiarity of the LF 
chiral transformation and explicitly give the null-plane 
chiral charge Q5 '. Quantum analysis, which is the main 
part of this paper, is developed in Sections IV and V. In 
Sec. IV, we demonstrate that perturbative and nonper- 
turbative treatments of the solution to the constraints 
give different description of the model. In Sec. V, we dis- 
cuss some physics consequences of the nonperturbative 
analysis. Especially, we resolve a problem of contradic- 
tion between the triviality of the null-plane charge and 
the nonzero condensate. We further calculate the masses 
of scalars in the symmetric and broken phases. Then 
we discuss the PCAC relations and nonconservation of 
the LF chiral charge Q\ Y in the broken phase. The last 
section is devoted to conclusion and discussions. 



II. THE MODEL 

Here we introduce the model (chiral Yukawa model) 
and summarize the standard knowledge on the chiral 
symmetry breaking in the conventional equal-time formu- 
lation. We also define our setup of the problem following 
the DLCQ method. 



A. Definition of the model 

The NJL model was first introduced as the simplest 
3 + 1 dimensional example which exhibits the dynamical 
chiral symmetry breaking In its original form there 
were two flavors, but a one flavor model 



C mh = ^ a {ip-m)m a + 



2N 



also breaks the chiral symmetry which exists in the mass- 
less case m = 0. We give an additional internal structure 
to the fermion independently of the flavor, and treat an 
iV-component spinor ty a in order to clarify the validity 
of approximation we use. From now on, summation over 
a = 1, . . . , N is always implied. Since this model is not 
renormalizable, we must specify a regularization scheme 
such as a cutoff to uniquely determine the model and to 
obtain finite results. 

In this paper, we discuss more general model with 
Yukawa interactions: 

C = V a {ip - m)* a + J^( d n< jd '*< J + d^d^ir) 



N 
2X 



{a 2 +7T 2 ) - (CT* a * Q + 7r* Q i7 5 * a ), (2.1) 



where a (71") is a scalar (pseudoscalar) boson with mass 
n/y/X and fj, is a dimensionless parameter. If one takes 
infinitely heavy mass limit for scalars fj, — ► 00, the dy- 



namical scalars become auxiliary fields a 



N 



7r = — -^■^ , Z75 , 3/ and the model goes back to the NJL 
model. 

In order to contrast with our LF calculation, let us 
briefly comment on the usual story of chiral symmetry 
breaking in the NJL model Q]. When m = 0, both of 
the lagrangian densities are invariant under the chiral 
transformation: 



cos 29 sin 29 
-sin 29 cos26> 



(2.2) 
(2.3) 



It should be reminded that this transformation is, of 
course, imposed on all the fields, which is, however, not 
the case in the LF formalism. This point will be dis- 
cussed later in more detail. The usual story is as fol- 
lows: The chiral symmetry breaks down spontaneously 
in a quantum level due to nonzero fermion condensate 
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(^>^>) 0. The most straightforward demonstration will 
be the mean field approximation with the concept of self- 
consistency. If one has TV-component fermion, we can 
justify the mean field approximation by the leading ap- 
proximation of l/N expansion. The self-consistency con- 
dition is a crucial key to the description of broken phase. 
This condition directly leads to the gap equation which 
determines the value of condensate and, equivalently, the 
physical fermion mass. As a result of symmetry break- 
ing, there emerges a Nambu-Goldstone boson. Since we 
do not have any fundamental scalar boson in the NJL 
model, the NG boson (pion) should be supplied dynami- 
cally as a bound state of a fermion and an anti-fermion. 
Mass of the pionic state indeed vanishes in the chiral 
limit. 

It will be helpful to co mme nt on the physics mean- 
ing of treating the model (2d). First of all, we should 



clearly distinguish our model from the linear a model of 
Gell-Mann and Levy ||l4| . Structurally our model resem- 
bles the a model in the sense that it consists of bosons 
and fermions interacting with each other via Yukawa cou- 
plings and has a continuous chiral symmetry. However, 
an important difference in our model is the absence of 
potential term for scalars. In the linear a model, it is 
the wine-bottle potential that induces the chiral symme- 
try breaking. Therefore the symmetry breaking occurs 
in the tree, level and the dynamical formation of NG bo- 
son can not be seen. What we obtain is only the sigma 
condensate (a) ^ 0. This naturally leads us to iden- 
tify the fermions with nucleons[j]. On the other hand, 
the fermions in our model should be regarded as quarks 
rather than nucleons. Indeed, as far as the leading or- 
der oi l/N expansion is concerned, the model shows the 
same behavior as the NJL model. For example, straight- 
forward calculations such as effective potential |L6| or 
mean-field approximation show that the chiral symme- 
try breaking (a) = — -4(# 1 ! r ) 7^ occurs in the one- loop 
quantum level (see Appendix A for more details). There- 
fore we do not consider our model as a special case of the 
linear a model and in order to remind this, we call it the 
"chiral Yukawa model" . 

It is also very important to view the NJL model as a 
low-energy effective theory of the chiral Yukawa model. 
The relation between two models is very similar to that 
between the Weinberg-Salam model and the Fermi theory 
of weak interaction: The chiral Yukawa model is renor- 
malizable and the fermions interact with each other by 
exchanging scalar or pseudo-scalar bosons. If we take 
the infinitely heavy mass limit for bosons /i/VA — > 00, 



1 However the fermions are sometimes treated as quarks. For 
example in the LF formalism, Carlitz et al. jl5| investigated 
the linear a model (i.e., with the wine-bottle potential) re- 
garding the fermions as quarks. But such treatment does not 
tell anything about dynamical chiral symmetry breaking and 
therefore should be clearly distinguished from our standpoint. 



then the theory reduces to the NJL model with non- 
renormalizable four-Fermi interactions. In this sense, 
the NJL model can be considered as a low energy ef- 
fective theory of the chiral Yukawa model. Such low en- 
ergy approximation will be valid when the momentum is 
much smaller than the mass of the exchanged particles 
p 2 <C M 2 /A- Lu this paper, we mainly treat the chiral 
Yukawa model for technical reasons, but what we even- 
tually want to know are results of the NJL model. There- 
fore even if we encounter divergences during calculation 
in the chiral Yukawa model, we only regularize them by 
some cutoff scheme and do not renormalize them. 

Finally, as far as we are discussing such 'low energy 
region', we do not have to worry about the problem of 
double counting of physics degrees of freedom. In the 
NJL model, the scalar and pseudo-scalar bosons are de- 
scribed as quark-antiquark bound states. On the other 
hand, we treat the scalars as physically independent de- 
grees of freedom in the chiral Yukawa model. So if we 
regard the fermions and scalars as quarks and mesons, 
there is the problem of double counting with which we are 
always confronted in treating chiral-quark type models 
p7| . However, in the low energy region, or equivalently 
for sufficiently large boson mass /1 — > 00, the scalars be- 
come "frozen" and do not behave as propagating degrees 
of freedom. 



B. Setup in the LF quantization 



We analyze the chiral Yukawa model (2d) in the 



DLCQ method and take special care of the longitudi- 
nal zero modes of scalars. In this method, we compactify 
the longitudinal space into a circle x~ £ [-L, L] with ap- 
propriate boundary conditions on fields. For scalars, we 
impose periodic boundary conditions at each LF time 



a(x = —L, x±) = a(x = L, x±) 
n(x~ = —L, x±) = Tt(x~ = L, xx) 



(2.4) 
(2.5) 



so that we can explicitly treat the longitudinal zero 
modes defined by 



a {x±) = ^- J dx <j(x), 
no(x±) = J dx~ n(x). 



(2.6) 
(2.7) 



Then the scalar fields are decomposed into the zero 
modes and the remaining oscillation modes: 



a(x) = a (x±) + tpa{x), 
ir(x) = 7T (:z; jj + ^(x). 



(2.8) 
(2.9) 



On the other hand, we impose an antiperiodic bound- 
ary condition for the fermion field, 



^ a (x- = —L, x x ) = -V a (x- = L, x x ). 



(2.10) 
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Here we must be careful about the boundary condition 
on the "bad component" of the fermion. As we discuss 
in the next section, if we decompose the fermion as 



(2.11) 



we find that i\>- ( "bad component" ) is a dependent field 
(see Appendix for the definition of A±). So the bound- 
ary condition on ?/>_ should be imposed consistently with 
the dynamics. For example, if we imposed the periodic 
boundary condition on ip_ and antiperiodic on ip + , the 
mass term tyty and the fermion's kinetic term became 
antiperiodic. This is not desirable as a term in the la- 
grangian and even not consistent with the scalar sector. 
Then how about the periodic boundary conditions for 
both of ip+ and ?/>_? In this case we have a dynamical 
zero mode of ?/;+, which is, however, not important to our 
problem because the chiral condensate will be related to 
the zero modes of a or ^vp. Periodic fermion will give 
unnecessary intricacy to the pr oblem . Therefore the an- 
tiperiodic boundary condition (2.1C) is appropriate. 



III. CLASSICAL ASPECT 

Classical analysis is necessary for specifying indepen- 
dent degrees of freedom. In this section, we determine 
the constraint structure of the model and define the LF 
chiral transformation. Chiral current and charge are ex- 
plicitly given. 



A. Constraints 

The system has three important constraints character- 
istic of the LF formalism: a constraint for 0_ (fcrmionic 
constraint) and two constraints for zero modes of bosons 
(zero-mode constraints). The Euler-Lagrange equation 
for itself is the fermionic constraint: 

id-^t = -(i7 X d± +m + a- ^75)7+?/'+ ■ (3.1) 

Also the zero-mode constraints for <7o(2;_l) and ttq(x±) 
are easily obtained from a; "-integration of the Euler- 
Lagrange equations for a and 7r, respectively: 



■--81 



l£_ 

N 



1 

«75 



0. 



where [ ]o denotes integration over x (see Appendix 
B). More explicitly, 







A 

N V2 



0] 



of 



-1 

«75 



,.«t (- 1 
\H5 



(3.2) 



These equations mean that gq and ttq should be repre- 
sented by other independent variables. If we take the 
fi — > 00 limit, the zero- mode constraints are reduced to 
zero-mode projected equations of the familiar relations 
a = and 7r = --^#£75*. Eventually the in- 

dependent degrees of freedom are nonzero modes of the 
scalars ip a , tp w , and the "good component" of the spinor 



The above constraints are, of course, derived from 



Dirac's procedure (see Appendix C). It is easily found 
that they belong to the second class. 



B. Chiral transformation on the LF 

Definition of chiral transformation on the LF is dif- 
ferent from the usual one Eqs. (2.2) and (|2.3[). This is 



because the identification of independent degrees of free- 
dom is not the same as usual. As we saw, <7o, ttq and ipt_ 
are dependent variables and should change as a result of 
transformation of the independent variables ip a , ip w and 
■0" . Therefore in the LF formulation, we impose the chi- 
ral transformation only on the dynamical variables: 



R{29) 



^7T 



(3.3) 



(3.4) 



where R(29) represents a rotation matrix defined in 
Eq. ( |2.3[ ). These are the definition of the "LF chiral 
transformation". If we find that tp- and (co,7ro) also 

as a result of (3.3) and (3.4), we can say that the "LF 



transform as 



3 ie ^-0_ and 



chiral tran sform ation " is substantially equivalent to the 
usual one (2.2) and (2.3). However, what is surprising 



about the "LF chiral transformation" is that the trans- 
formation ( |3.3| ) is an exact symmetry even for massive 
fermion as far as interaction is absent Jig] ]. So it will be 
interesting to check whether the "LF chiral transforma- 
tion" in our model is exact or not when a mass term is 
present. 

In order to see the transformation property of the de- 
pendent variables, let us solve the constraints classically. 
This means that we completely ignore the ordering of the 
variables which becomes a burdensome but important is- 
sue in a quantum treatment. The fermionic constraint 
(|3.l| ) which was originally a complicated relation in the 
purely fermionic NJL modelQ, is now easily solved owing 
to introduction of scalars. The zero-mode constraints 



2 It is difficult but possible to solve the fermionic constraint 
in classical treatment where we just treat the spinors as Grass- 
mannian numbers. The exact solution obtained is highly non- 
local and complicated 12 1. 
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are also solved formally. Explicit form of the solutions is 
given in Appendix C. Now we find the transformation of 
the zero modes and subsequently that of i/ 1 - ■ 



1. Massless fermion 

Let us first consider the massless fermion case. When 
m = 0, it is easy to see that the transformation (3.3) and 



(3.4) induces the following: 



R{26) 



(3.5) 



(3.6) 



This is identical with the usual chiral transformation. 
Therefore it i s sh own that when m = the fields do 
transform as ( |2 . 2| ) and (|]^) even on the LF at the clas- 
sical level. 

Now that we know all the transformation laws, it is 
straightforward to construct the Noether current and 
charge. The lagrangian with m = is invariant under 
the LF chiral transformations. Form of the LF chiral 
current is equivalent to the usual one, 



(3.7) 



However, do, 7i"o and ?p- in (3.7) should be understood as 
solutions of the constraints. On the other hand, the LF 
chiral charge Q\ F — J dx~ J d 2 x±_ jg" does not include 
the constrained variables: 

Q\ F = Jd 3 x -V2^W+ + ^(^d-^ a - (p a d-ip T 

(3.8) 



which is consistent with the fact that Q^ F is a generator 
of the chiral rotation for independent variables. Trans- 
formation of other dependent fields should be obtained 
through the change of dynamical variables. 



2. Massive fermion 

The massive fermion case is much more complicated. 
As mentioned before, an astonishing fact of the "LF chi- 
ral transformation" is that it is an exact symmetry even 
for a massive free fermion |l|. When the mass term 
is present, the "bad com ponent" of free fermion does 
not transform as Eq. ( |3.6| ). Subsequently, the associated 
Noether current has an extra term proportional to the 
bare mass 



75 t5 



im(ip + + ip_)j fi j 5 —j + ip + . 



(3.9) 



Nevertheless, the divergence of the current turns out to 
be zero = due to the cancellation between the 

first term {d^J^ = — 2im^75^) and the second term. 
This should be compared with the usual current in the 
equal-time quantization: = — ^'7^75^', d^ J^ — 

— 2ir7i5'75^', which also holds for interacting theories and 
is intimately connected with PCAC relation. Note that 
the LF chiral charge Q\ F = J d 3 xJ§ is equivalent to that 
in the massless case due to = . This is natural be- 
cause the LF chiral transformation is defined irrespective 
of the mass term. 

Now, how about the chiral Yukawa model? Using the 
solution for m ^ (see Appendix C), the infinitesimal 
chiral transformation of <7q, ttq and are given as fol- 
lows: 



f Scr 
\5ir 



26tt 
-26a 



( ° 



2z0 75 m(l + C)i^-7 + ^+, 
2 i d- 



(3.10) 
(3.11) 



where 



V2N V ±y 



and T>(x±) is the transverse differential operator defined 
in Appendix C. It is evide nt t hat th e d ependent fields 
do not transform as Eqs. ( j3~5| ) and ( |3.6|) . As a result, 
the Noether current ( |3.7j ) also gets modified by the term 
proportional to m. The explicit form of the current j§ is 
very complicated but we can see that the -I- component is 
equivalent to that of ( 3.7). Therefore the LF chiral charge 
is given by Eq. (3.8) even for massive case. However, 
contrary to the massive free fermion, the divergence of 
the LF current <9 M j^ does not vanish due to nontrivial 
interactions. 

One of the lessons suggested by these observations 
is that when we investigate physics related to massive 
fermions (e.g. PCAC relation), we had better treat the 
current defined by the massless fermions rather than 

j?. This is clear for the free case: The true LF chiral cur- 
rent ( |3.9| ) defined for the massive fermion vanishes if 
we take the divergence while that for massless fermion 
gives the usual relation. This is true of the chiral Yukawa 
model. The divergence of the current (3.7) is given as 



(3.12) 



while <9 M j^ is very complicated. It will be very difficult 

to discuss the PCAC relations etc. by using j'jj. It is not 

quite clear if analysis of makes sense. Instead, even 
for the massive case, we treat the current ( |3.7D to discuss 
the physics such as PCAC. This point will be discussed 
later again in Section V. 
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IV. QUANTUM ASPECTS 

In the classical analysis, we formally solved the con- 
straints in order to find the LF chiral current and charge. 
When m = 0, the resulting hamiltonian is chiral sym- 
metric and we do not have any symmetry breaking term. 
Therefore even if we go to quantum theory with such 
hamiltonian, we will not be able to describe the chiral 
symmetry breaking. Certainly it might be possible that 
we could find a broken phase hamiltonian by adjusting 
the operator ordering, but such procedure seems unnat- 
ural and tricky. Instead, we quantize the model before 
solving the constraints. This means that we perform the 
Dirac quantization for constrained systems. After that, 
the constraints are solved quantum mechanically with a 
care of the operator ordering. The same route has been 
traced by many people who tried to describe the sponta- 
neous symmetry breaking of simple scalar systems |5| . 

Calculation of the Dirac brackets in our system is 
a very complicated task. However, the Dirac brackets 
between dynamical variables turn out to be standard 
ones: Quantization conditions for the dynamical vari- 
ables lpo-tPtv, and ip+ are 



. m 2 i 



IT) 



S (2) 



(4.1) 

y±), 



= -±=(A + ) afj S ab 5(x- - y-)8&>(x ± - y±), (4.2) 

where £ and rj stand for a or it, and a, (3 — 1, • • • , 4 are 
the spinor indices. The sign function e(x~) is defined 
in Appendix B. The other commutators between dynam- 
ical variables are zero. Note that these conditions are 
irrespective of the phase of the model because they are 
independent of the interaction. 

Mode expansion of the fields at x + = reads 



V>+(x) 



d 2 pj_ 
■ , (2tt) 2 

x {a v (Pn,P±) e~' ipx + a^(p+,p 

= _ J_ \p f d 2 p 

2L ^ 



e ipx }, (4.3) 



1 1 

" 2 ' 2 



Y, {w(h)b a (p+,p ± ,h)c- 



h=±h 



+ w(-h)d a \p+,p ± ,h) e J P x }, (4.4) 

where px = Pn%~ — P±x^ and p+ = nn/L. The spinors 
w(±h) depend only on the helicity h It follows that 



= {27r) 2 (2L)2p+S nm S{ P± - q ± )S (v , 
{b a (pt, P± ,h), V)} 

= 2p+|— 8 nm S(px - q±)5hh>, 
Ztt 

{d a (p+,px,h), d b \q+,q ± ,h')} 

= 2p+|— 8 nm S(px - q±)5hh>- 
Ztt 

It is important to note that both of the above mode 
expansions are independent of the mass. (The spinors 
w(±h) are independent of mass. This is clearly shown in 
the Appendix of Ref. §.) This means that if we calcu- 
late n-point Green functions at fixed time x + = 0, they 
will become independent of the value of mass, which is 
not a correct result in general. This undesirable situation 
is known as one of the pathological properties of the LF 
formalism which needs great care for obtaining correct 
results p8| . Indeed, as we will see later, to remedy this 
problem is indispensable to get a meaningful gap equa- 
tion. In many cases, loss of mass information is cured by 
a carefully chosen infrared rcgularization. It should also 
be commented that the mass-information loss is certainly 
not a desirable feature, but we will find its usefulness in 
describing the broken phase physics. Anyway, we must 
pay great attention to the fact that naive mode expan- 
sion of the (scalar and fermion) fields is independent of 
the value of mass. 

The Fock vacuum |0) is defined by 



a<r(p£) PJ-)|0) = a n (p+, p±)\0) 

= b a ( P +,p ±: h)\0)=d a (p+ 1 p ± ,h)\0)=0, 



(4.5) 



for n > 0. It is worth while emphasizing that the vacuum 
of this system is really the Fock vacuum since we have no 
dynamical zero modes. Because of the p + conservation, 
the normal-ordered chiral charge Qg F always annihilates 
the vacuum: 



Q\ ¥ \0) = 0. 



(4.6) 



It has been known that any light-like charge Q LF au- 
tomatically leaves the vacuum invariant Q LF |0) = 
whether or not it generates a symmetry [ fl9|| . 

In a quantum theory, operator ordering becomes an 
issue. Let us comment on the problem of operator or- 
dering and clarify our stance toward it. Since the (zero- 
mode and fermionic) constraint equations are generally 
nonlinear relations among operators, their solutions de- 
pend on operator ordering. We must select an appropri- 
ate operator ordering. Then, what can be the criterion 
for this problem? In many papers discussing the spon- 
taneous symmetry breaking in DLCQ, the Weyl order- 
ing is adopted on general grounds. However, it is not 
clear whether the Weyl ordering in constraint equations 
makes sense because they include both independent and 
dependent variables. The most reliable criterion for de- 
termining the operator ordering will be as follows. Be- 
fore solving the constraint equations, we can calculate 



G 



the Dirac brackets between independent variables and 
dependent ones (e.g. [ao,ip+] — ■ ■■), which are terribly 
complicated in our model and we do not display them 
in this paper. Here we already have to specify the op- 
erator ordering. On the other hand, we can solve the 
constraint relations with the above ordering and obtain 
their solutions such as oo = cro(ip a , <Pn, ip+)- Now we 
can calculate again the commutators between the solu- 
tions (i.e. dependent variables) and independent vari- 
ables (e.g. [po(y?o-, y^ , ip+] ) using simple commu- 
tators Eqs. (4.1) and Q4.2| ). The results should be iden- 
tical with those of the Dirac bracket. In other words, 
we must find out such operator ordering that will give a 
consistent result in the above sense. This should be the 
criterion for an appropriate operator ordering. However, 
as you expect, to find such ordering in our model is an 
extremely difficult task. So practically we just work with 
several particular orderings and compare the results. In 
our actual calculations, we treat two specific orderings 
and check whether the results depend on the ordering 
or not. To find a consistent operator ordering should be 
examined in much simpler models. 

In the following, we will solve the zero-mode con- 
straints in two different ways: perturbative and non- 
perturbative methods. To solve the constraint is signif- 
icant to describe the symmetry breaking on the LF. To 
see this, let us decompose the longitudinal zero modes 
into c-number parts and normal-ordered operator parts, 



(c) 

Co = cr + <t, 

(c) 



(OP) 
> 

.(op) 




(4.7) 
(4.8) 



If the c-number part of the solution is nonvanishing, it di- 
rectly means nonzero condensate: (0 |<x| 0) = Oq 7^ an d 
(0 1 7T 1 0) = 7Tq C ' ) 7^ 0. Therefore to find such nontrivial solu- 
tion is necessary to describe the symmetry breaking. We 
explicitly demonstrate that perturbative solutions can- 
not lead to chiral symmetry breaking while nonpertur- 
bative solutions give nonzero vacuum expectation value 
for a. In both cases, the fermionic constraint is formally 
solved (as in Eq. (|C2D) and inserted into the zero-mode 
constraints. 



■»-£(£)>■ 

n— v / 



(4.10) 
(4.11) 



and the dynamical variables are treated as 0((-£-)°). 
Inserting the above expansions into the constraints 
Eqs. (3.1) and ( |3.2] ) with the natural ordering and com- 
paring the same order of perturbation, we obtain the so- 
lution order by order. For example, the lowest order 
solutions are 



CT (0) _ (0) _ 



(4.12) 



o(0) 



1 1 

2ibT 



(ij x d±+m + ip (7 -wyfe)7 + ^+- (4.13) 



Higher order solutions are given in Appendix D. 

The chiral transformation of the perturbative solution 
in the massless case can be inductively checked. First, 
it is easy to see that the 0-th and 1st order solutions ro- 
tate symmetrically under the LF chiral transformation 
( |3.3D and (3.4). If we suppose the n-th order solutions 
rotate chirally, the n + 1-th order solutions also behave 
the same. Eventually the perturbative solution trans- 
forms symmetrically under the LF chiral rotation ( |3.3| ) 
and (O), 



(4.14) 
(4.15) 



[QJV] =-2wr, [Q™,n]=2ia 



Of course the hamiltonian [cf) Eq. (CI)] is also invariant 



[Q^ F , H] = which is the same as the classical analysis 
in the previous section. 

What is most important is that the vacuum expectation 
values of the perturbative solutions vanish in all order 
of perturbation (0|<7o|0) = (0|7To|0) = 0. This is easily 
verified by using (0|v? CTi7r |0) = and (O|^_(9±<0+|O) = 0. 
Therefore we are in a chiral symmetric phase: 

(OjcrjO) = <0|tt|0) = 0. (4.16) 



A. Perturbative solutions to the zero-mode 
constraints 

Let us solve the zero-mode constraints using pertur- 
bation in terms of the coupling constant A. Since A is 
a dimensionful parameter, we introduce some scale A cr 
which is much larger than A (A cr 3> A). We regard A cr as 
a critical coupling of the symmetry breaking which will 
be determined later. Now we expand the constrained 
variables as follows 

n— x ' 



B. Nonperturbative solutions to the zero-mode 
constraints 

We next solve the zero-mode constraints using the 
mean-field approximation. In the following, we work 
with a particular operator ordering though the result is 
the same as others as far as we discuss the leading or- 
der of 1/N expansion. The following ordering greatly 
reduces our calculation. Substituting the solution of the 
fermionic constraint into the zero-mode constraints and 
rearrange the ordering, we obtain 
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1 fl 2 1 

+ W 2 ~N 2~L 
x {^(s) 



L 

L 
1 



-L 



2/ 
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7^V+(y) 



-1 



m + a{y) 

-n{y) 

m + a(y) 
h.c, 



(4.17) 



where ip+(y) = y~ and similarly for cr(y) and 

it{y) The operator ordering here is different from that in 
the perturbative treatment. However, one can show that 
the previous perturbative result does not change with 
the above ordering. That is, the perturbative solution 
with the above ordering does not lead to chiral symmetry 
breaking. 

Let us first determi ne t he c- num ber part of the zero 
modes defined by Eqs . (4.7) and (4.g). We saw in the clas- 
sical analysis that (oo, ttq) rotates chirally in the massless 
case. Therefore we choose 



0. 



(4.18) 
(4.19) 



Taking a vacuum expectation value of the zero-mode con- 
straint for a greatly simplifies the calculation, which is 
an advantage of our specific choice of ordering: 



£_ (c) 



1 V f , ( c h 



Introducing M defined by 
M = m 4 



1 

id- 



7 (c) 



(4.20) 



(4.21) 



and eval uating the vacuum expectation value in 
Eq. ( 4.20 ) by using the mode expansion, we have 



M 



1 A 



;M< 



V2N 



n=l/2, 



Ap 4 
pt 



(4.22) 



where Ap + = tt/L. The operator form of the right-hand- 
side suggests us to identify as 



M I 

7=2 \° 
= <o|[*Sf*s f ] |o> > 



T id- 



(4.23) 



where ^>\ { is a fermion with mass M, 



^ = i-L(i r L J _ + M) 7 +^ . 



2 id- 



(4.24) 



(V'-a/ is the "bad" component of \&5Lf.) Therefore it is 
natural to consider M to be the physica l ferm ion mass. In 
other words, the identification in Eq. (4.23) corresponds 
to the self-c onsistency condition. 

Eq. (4.22) should be the gap equation by which we 



(c) 

can determine a and equivalently, the physical fermion 
mass M. However, it is not evident to regard it as the 
gap equation because Eq. (4.22) in the chiral limit m — > 
cannot give nonzero M. The same situation was observed 
in our previous work on the Gross- Neveu model [p"l| . 
As was discussed in Ref. ]Tl| , if we want a meaningful 
gap equat ion, we must supply mass information so that 
Eq. possesses a nontrivial solution M ^ in the 

chiral limit when we regularize the divergent summation 
over n. The need of the mass dependence in Eq. ( |4 . 2 2|) i s 
readily understood from the identification in Eq. ( ^.23[ ). 
Indeed, one can easily check that (0|*&m ^Af |0) /M should 
explicitly depend on M in the equal-time formulation. 
This is a typical example of the "mass-information loss" 
on the LF |18[] which must be repaired properly for ob- 
taining correct results. 

It may be p ossibl e to regularize the divergent sum- 
mation in Eq. ( 4.22 ) with, say, a heat-kernel damping 
factor jyj, but such calculation is complicated and not 
tractable. Instead, we introduce some cutoff that ren- 
ders the divergent summation into finite one. Such cutoff 
should be introduced so that the result correctly depends 
on the mass M. Here for simplicity, we adopt a cutoff 
which eventually reduces to the parity invariant (PI) cut- 
off < A H . From the dispersion relation and the PI 
cutoff, we find that the momentum region is restricted to 
(M 2 + pi)/2A <pt = nn/L < A. Therefore we set 



niB, < n < n vv , 



(4.25) 



where nm and nuv are nearest half- integers to ^ — ^ p± 
and —A respectively. If we use the approximation for a 



large half-integer n, J2n=i/2 
summation is approximated as 

nuv-l 



Inn + ln4e 7 , the 



E 

n=l/2, 



pi 



i^ln^X 



In 



2A 2 



n=niR+l 



M 2 + p 2 ± ' 



(4.26) 



for fixed A, M and and sufficiently large L. Of course 
there is "finite volume effect" for finite L, but we finally 
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take the infinite volume limit and the finite volume ef- 
fect is expected to be small as far as L is large enough^. 
Eventually Eq. (4.22) becomes dependent on the mass M 
and can be considered to be a gap equation 



A 2 



M 



M - m = \M—{2-—[l + \n—)\. (4.27) 



A 



2A 2 



M 2 



Form of this gap equation is different from those with 
familiar cutoff schemes such as the three or four momen- 
tum cutoff H , but our gap equation behaves exactly the 
same as usual. Indeed, even in the chiral limit m — > 0, 
this equation is a nonlinear equation for M and when 
the coupling constant A is larger than the critical value 
A cr = 27r 2 /A 2 , there is a nontrivial solution M = Mq ^ 
(see FIG. 1). This also means that the zero mode of 



Further- 



(c) 

a has been determined as <t = M 
more, one should note that the gap equation (4.27) is 
independent of the value of So we can regard the fi- 
nite [i result (a) = M — m as the result for infinite /x; 

= M—m. (Remember that the /z — > co limit of 



' N ' 



(3.2) is Co = — ^[^^]o«) Therefore the chiral symmetry 



breaking occurs for arbitrary value of /x in the mean-field 
approximation. This is consistent with the result of the 
conventional equal-time quantization (see Appendix A). 

On the other hand, there is a trivial solution M = 
(when m = 0) even for A > A cr and if we select this 
solution the resulting theory becomes chiral symmetric. 
Then there comes a problem which solution should be 
physically realized. Unfortunately, comparison of the 
vacuum energy for both phases does not tell anything 
about this problem because the vacuum energies turn 
out to be the same. If we found the consistent operator 
ordering as discussed before, we could estimate difference 
of the vacuum energies and determine the physically re- 
alized phase. Even without such calculation, however, 
we can say that the symmetric solution is excluded for 
A > A cr . This is because there emerge tachyonic modes 
and the system becomes unstable if we select a trivial 
solution for A > A cr . This will be again discussed in Sec- 
tion V-B. So we deal with only the nontrivial solution for 
A > A cr and do not consider the symmetric solution. 



3 Finite volume physics itself is intriguing. For example, sim- 
ilarly to the equal-time calculation fed] , if we make the volume 
smaller and smaller, we will meet a critical length L CT beyond 
which the chiral symmetry never breaks down. Moreover, if 
we could determine the L dependence of the physical mass 
M(L), it would serve as a prediction for the limiting behavior 
of eigenvalues in the numerical DLCQ calculation. Neverthe- 
less, such finite volume physics is outside the scope of this 
paper and we do not discuss it anymore. We always assume 
L sufficiently large and ignore the finite volume effects. 




0.5 



2.5 



FIG. 1. Fermion's dynamical ma ss M p in the chiral limit 
as the solution of the gap equation (4.27). There is a nonzero 
solution for A > A cr . 



Comments on other cutoff schemes are in order. We 
find a nontrivial equation for M by using the PI cutoff. 
It was crucial to include the mass information as the reg- 
ularization. However, we have to be careful in setting 
the cutoff. Any cutoff scheme which holds mass does not 
necessarily lead to a physically sensible result For 
example, a two dimensional PI cutoff M 2 /2A < p + < A 
with a transverse cutoff |pj_| < A gives a wrong re- 
sult. The resulting gap equation erroneously predicts 
that there is no symmetric phase. It seems important to 
introduce a cutoff with some symmetry considerations. 
Indeed, the three momentum cutoff respecting the three 
dimensional rotation and the PI cutoff || in our case 
predict the existence of the critical coupling constant. 

Now let us determine the operator parts of the zero 
modes by the mean-field approxi matio n. We approxi- 
mate the nonlinear terms in Eq. (4.47) by using AB ss 
A(B) + {A)B — (A)(B), where the expectation values are 
taken with respect to the Fock vacuum. We further ne- 
glect contribution from the oscillating modes of scalars 
<y9<r = ^P-k = 0. Then the operator parts are given by 



r (°p) 
>p) 



• ( m ZM ~ d l) 



\5l a 
* M 



1 

«75 



where a c-number quantity to| m is defined as 



'■ZM 



/£ 
A 



J o 
(4.28) 



(4.29) 



N v/2 V 



The numerical value of 



l ZM 



is calculated if wc utilize 



the gap equation: m| M — ^jp- Inserting the c- number 
and operator parts of ao and 7To into the solution of the 
fcrmionic constraint, we have 
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(op) 




75 7 + ^, (4-30) 



where tp-M is given in Eq. ( 4.24| ). 

To understand what we did above, let us consider the 
relation between our operator ordering and the l/N ex- 
pansion. We have obtained an equation for the c-number 
part of ag (the gap equation) just by taking the vac- 
uum expectation value of the zero-mode constraint even 
without recourse to the l/N expansion. This simplic- 
ity in obtaining the gap equation is mainly due to our 
specific choice of the operator ordering. As has been 
commented before, if we take other orderings, our calcu- 
lation becomes terrible because of complicated structure 
of the Dirac brackets between constrained zero modes 
and physical variables. However, as far as the leading 
term of the l/N expansion is concerned, the commuta- 
tor [<7o,'0+] turns out to be of the order of 0(1/N) and 
we can ignore the effect of ordering. Furthermore, the 
approximation neglecting the scalar oscillating modes is 
also justified b y th e l/N expansion. From the quantiza- 
tion condition (4.1), we find ip a is 0(N~ 1 / 2 ) whereas Co is 
O(N ). These considerations justify that our mean-field 
calculation with the specific operator ordering is correct 
up to the leading contribution of the I/TV expansion. 

Before ending this section, it will be better to point 
out the "merit" of the mass-information loss. Certainly 
it was a demerit in deriving the gap equation, but this 
property gives a very important benefit to our frame- 
work. The fact that the mode expansion is independent 
of the value of mass in turn means that the Fock vacuum 



defined by (4.5) keeps invariant even if we change the 
value of mass. We do not have to perform the Bogoli- 
ubov transformation on the vacuum depending on the 
change of mass. Therefore the LF vacuum is invariant 
even after the fermion acquires dynamical mass M ^ 0. 



V. PHYSICS IN NONPERTURBATIVE REGION 



A. Null-plane chiral charge vs chiral condensate 

In the equal-time quantization, the broken vacuum 
does not possess the chiral symmetry Qf T |0) ^ 0. 
The Nambu-Goldstone phase is characterized by a 
nonzero condensate of the order parameter ('J"!') and 
the more strict expression of Qf T |0) ^ is a rela- 
tion (0|[Qf T ,*i 75 *]|0) = (0|J d 3 x \j°(x),m 75 *}\0) = 
2z (0 1 \E r \T/ 1 0) 7^ 0. Therefore there is no inconsistency 
between these two relations. On the other hand, re- 
member that the light-like charge always annihilates 
the vacuum. This implies that if a similar relation 

[Q\ %*i75*] =_2i** held on the light front in the 
broken phase (0|'J n I'|0) ^ 0, it would immediately con- 
flict with the triviality of the chiral charge Q\ F \H) = 0. 
In the following we resolve this seemingly inconsistent 
situation. 

The chiral charge Q\ F defined by Eq. ( |3.8| ) annihilates 
the vacuum and generates the chiral transformation for 
the independent variables irrespective of the symmetry. 
Indeed, we find 



[Of ,</>+] = 75 V + - 



2iip a 



(5.1) 
(5.2) 



These are the fundamental laws of the chiral transforma- 
tion. Any transformation of the dependent fields cto, ttq, 
and tp- should be derived from them. 

In the broken phase A > A cr , the gap equation has a 
nontrivial solution M ^ and the fermion behaves as a 
massive fermion with the dynamical mass M. First of 
all, let us view the chiral transformati on o f the massive 
fermion operator ^>m defined by Eq. (4.24). The result 
is already unfamiliar to us: 



M- 



1 1 



M> a M = -2M l5 -—^% 
2 id- 



(5.3) 
(5.4) 



In this section, we discuss some physics consequences 
of our method. Firstly, we explicitly demonstrate how 
the triviality of the null-plane chiral charge Q\ F and the 
nonzero chiral condensate reconcile with each other. Sec- 
ondly, masses of the scalar and pseudoscalar bosons are 
calculated from the lagrangian for both phases. Finally, 
we d erive the PC AC relation for the chiral current j£ 
[Eq. (|3.7[)1 and discuss the nonconservation of Q\ ¥ . 



The second term A^j^ does not exist in the equal-time 
quantization. Only if M ^ 0, this is equivalent to 
the usual transformation. Using this result, the chiral 



(°P) 



tran sform ation of a h 
Eq. ((Oof) are given as follows: 



in Eq. J4.2SQ and ip_ 



Q 



LF (op) 



Q7, 



7T, 





(op) _ 




-2in, 



(op) 



Act, 



(op) 



2ia, 



(op) 



Att 



(op) 



[or, r-] =^- + ( 1 +2M^o 



(o P ; 



A* Q M , 



(5.5) 
(5.6) 

(5.7) 



where 



4 If one takes fj, — > oo, one will be convinced that [cro, = 

[-±m}o,r+]=o(N- 1 ). 



Act, 
Att, 



(op) 

o 

(op) 



2i 



N m| M 



M 

71 



-i ( *75 



id. 
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-1 



*m|0) 



-tpl +h.c. 
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One can easily show that Act, 



(°P) 



is zero due to 



(0\^ M~i J 5^ m\0) — and antisymmetry of the sign func- 
tion e(x~ — y~). Therefore the terms involving A^ M 

and Att^ 9 ^ are the extra compared with the usual chi- 
ral transformation. They do not vanish even in the chiral 
limit. These extra terms are direct consequences of being 
dependent variables and the dynamical generation of the 
fermion mass M. Unlike the equal-time calculation, the 
chiral transformation of the full field variables becomes 
model- dependent in general because a part of the vari- 
ables are constrained and the information of interaction 
inevitably enters the transformation law of constrained 
variables through the solutions. 

Due to the modification of the transformation law, we 
can avoid the inconsistency. The transformation of the 
full 7T field is given by [Q% F , tt] = 2i(a [ ° p) +^)+a4° p) ^ 
2ia. The vacuum expectation value of this equation gives 
a consistent result: 



(0|[Qf , tt]|0) = 2i(0 



>p) 



Att, 



(°P) 



= 0. 



(5.8) 



Now it is easy to obtain the transformation of \N75\E'. 
Our final result is 



(5.9) 



^«t2_ (*a4 op) + 2M 



1 1 
2~id~- 



7 + ^° + h.c. 



In addition to the first term that is equivalent to the usual 
result, we have nonvanishing extra terms. However, if one 
takes the vacuum expectation value of this equation, such 
extra terms should exactly cancel the first term 2i (^>^) ^ 
0. It is indeed the case and the explicit evaluation of the 
r.h.s. gives a consistent result 



<0|[Q^ F ,§z 75 *]|0) =0. 



(5.10) 



Here we neglect the term [Att^° p) , ip^} ~ 0(1/N). Thus 
we checked the consistency between the null-plane charge 
Q\ F and the chiral condensate (a) = — (<i n I') ^ upto 
the mean-field level. 

As the result of these unusual chiral transformation, 
even the hamiltonian loses the chiral symmetry in the 
broken phase: The commutator [Q$ F , H] ^ is directly 
evaluated exactly in the same way as above. This means 
that the LF chiral charge is not conserved in the bro- 
ken phase. It should be emphasized that the violation 
is proportional to the fcrmion's dynamical mass M and 
thus does not vanish even in the chiral limit. It is very 
interesting that the chiral symmetry breaking in the LF 
formulation is expressed as an explicit breaking. The im- 
portant difference, however, is that usual explicit break- 
ing does not accompany the gap equation, while in our 
case the gap equation plays a very important role in many 
aspects. The nonconservation of Q\ ¥ on the LF in the 



broken phase has been discussed by several people in re- 
lation to PCAC []l5 22 2^]. Particularly, a similar situa- 
tion to our conclusion (i.e. the non-conservation of the 
null-plane charge in the DLCQ method) was found in 
the broken phase of the scalar model |2j| . The problem 
of nonconserving charge should be intimately connected 
with the divergence of the chiral current. In Sec. V-C, we 
will again meet the non-conservation of Q\ F as a result 
of the PCAC relation and peculiar behavior of the pion 
zero mode in the chiral limit. 

Let us turn to the symmetric phase where the cou- 
pling is large but slightly less than the critical value 
A ~ A cr . In this region, we use the symmetric solution 
of the gap equation. If we restrict ourselves to the chiral 
limit m = 0, the solution is just a trivial one M = 0. 
Transformation law in this phase is obtained by simply 
substituting M = into the above results. Therefore in 
the leading order of 1/N expansion, all the dependent 
fields transform in a chiral symmetric way: 



(i 



LF r _ 



(op) 



[Q 



= — 2in, 



(op) 



= 2ia, 



(op) 



(5.11) 
(5.12) 



With these commutators, we find also [Q^if] = 0. 
Thus the LF chiral charge is conserved in the symmetric 
phase as expected. 



B. Masses of the scalar and pseudoscalar bosons 

Here we calculate the masses of scalars which are 
bound states of fermion and antifermion. Unlike the NJL 
model, we have 'dynamical' scalars in the lagrangian. 
Therefore it is convenient to evaluate the "pole mass" of 
the scalars directly from their propagators without con- 
sidering bound-state equations. The procedure of calcu- 
lating the pole masses is as follows. First we insert the 
broken or unbroken solutions ao, ttq, and ip_ into the 
original lagrangian. For simplicity, we ignore the finite L 
effect. This is because we are only interested in the ef- 

(c) 

fects of the condensate Cg and the nonzero constituent 
mass M. Next, reading the fermion propagator from the 
lagrangian, we calculate the scalar propagators A 7r (T (A; 2 ) 
upto one loop of the fermion. Finally, the pole masses 
are obtained from the equation A~\(k 2 = m\ n ) = 0. 



1. Broken phase 

Let us first consider the broken phase A > A cr . Insert- 
ing the broken solutions into the lagrangian, we have 



-V2^id+^ a + 



n 



x {V>-m + g JdZ ( & ~ ^7s)7 + V>+ } 



N 

v 



(c^a + S^tt) - ^ {(a< c) + a) 2 + tt 2 } 



M 

(a + i7r 75 ) 7 +^-{(a + i7r 75 )^} , (5.13) 

where we used the notation a(x) = + <r(x) and 
7r(x) = 7r(x), and is defined by Eq. ( 4.24 ). Instead 
of the fermion propagator for + component S++(p) = 
\/2k+p- / (p 2 — M 2 + it), it is convenient for practical 
calculation to define the propagator for 



m- * +M 



P 2 -M 2 



ic 



(5.14) 




FIG. 2. Fermion's one loop contribution coming from the 
Yukawa interaction. The solid line is for the fermion ^ m, the 
dashed line for 7T or <j. 



where p M = (p + ,p~ = (p 2 + M 2 )/2p + ,p±) is the on- 
shell four momentum [M. Note that this partially on- 
shell propagator S(p) is different from the usual fermion 
propagator S(p) by an instantaneous part S(p) = S(p) — 

^qr, which arises from the bad component ip_ as the 
solution of the fermionic constraint. 

Scalar and pseudoscalar propagators A CT)7r (fc) with 
fermion's one loop quantum correction are given by 



N I A" 1 (A;) 



k 2 1 

— --+ F Yu kawa(fc) + F inst (k), (5.15) 

fi A 



where 

-FVukawa(fc) = I [d 4 p] tr 



l l5 \ - 



S( P ) 



(5.16) 



comes from the Yukawa interaction (FIG. 2) and 



F inst (k) = -- I [d 4 p] tr 



Sip) 
lis \ 1 



lib 
1 



p+ - k+ 



«75 



H5 

1 



1 J p+ + k+ \ 1 



(5.17) 



from the instantaneous interaction (FIG. 3). The inte- 
gration measure is given by 



[A] 



d 2 Pl _ 
i{2ir) 2 



dp+ 
~2w 



dp 



Summation over the longitudinal discrete momenta p+ is 
approximated by integration. 




FIG. 3. Fermion's one loop contribution coming from the 
instantaneous interaction which is represented by the vertical 
solid line. 



For simplicity, we put k±_ — 0. Using a parameter 
x = p + /k + , the propagators A W (T (fc) are expressed as 



N\A-i(k) ..2 



rn 



T{k 2 ,M 



A* M 

1 r<fi 

dx 



d 2 p_ 



ZXI <-T(k 2 ,M 2 ) 
1 



k 2 

k 2 - 4M 2 



(2tt) 3 p\ + M 2 - k 2 x(l - x) 



(5.18) 



where we have utilized the gap equation (4.22) and 



m 



ZM 



fjrm/XM. The physical masses are determined 
from the equations A~\(k 2 — m 2 ^) = 0. Since the in- 
tegral in T{\^,M 2 ) diverges, we must specify a cutoff. 
Here we use the "extended PI cutoff" 11(1: 



(5.19) 



where i denotes the particles of the internal lines. This 
cutoff is a natural generalization of the naive PI cutoff 
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p± < A and can be applied to multiple internal lines. In 
our case, the extended PI cutoff becomes 



Pi 



M 2 



Pa 



M 2 



l-x 



< 2A 2 



The explicit form of .F(fc 2 , M 2 ) with this cutoff is given m 
Appendix E. Then we obtain highly complicated nonlin- 
ear equations for m 2 a which are also shown in Appendix 
E. First of all, it is almost trivial that the equation for 
m n has a solution m, = in the chiral limit. For small 
bare mass m ~ 0, we find 



XMn 



1 



^"EPl(fc 2 -> 0, Mq) 



0(m 2 ), 



(5.20) 



where Mq is a fermion condensate in the chiral limit^] and 
^EPi(k 2 — > 0, Mq) is independent of m (see Appendix E). 
Therefore we have checked that the mass of the pion goes 
to zero in the chiral limit and now we can identify it with 
the NG boson. 

The mass of a is determined in the same way. For 
example, in the chiral and heavy mass limit m —* 0, [i — > 
oo, one can easily find a solution m 2 — (2M) 2 which is 
known to exist in the NJL model in the chiral limit. 

Our result of the pion mass ( fT2C| ) satisfies the Gell- 
Mann-Oakes-Renner (GOR) relation when /i — > oo: 



m 2 / 2 = -4m(*f) 



(5.21) 



where f v is the pion decay constant. To see this, let us 
calculate /„• explicitly. Rewriting the pion's propagator 
in the chiral and heavy scalar limit as 



A w (fc) = 



k 2 — m 2 
Z^ = [NT EP1 (k 2 



0,M 2 )]- 



(5.22) 
(5.23) 



therefore 

U = 2M [NT EPI (k 2 -> 0, M 2 )} 1/2 = ^ . (5.25) 

v Zn 

Using this result and m 2 = jjj-NZk, we finally confirm 



the GOR relation fl5.2l| ) 



2. Symmetric phase 

Next let us consider the symmetric case. Using the 
symmetric solution to the zero-mode constraint, we can 
evaluate the masses of a and 7r in the symmetric phase. 
In the chiral limit, we have 



'ZM 



N 



where the zero-mode mass 



T(k 2 ,M 2 = 0)k 2 (5.26) 



'ZM 



in the symmetric phase 
is expressed differently from that in the broken phase: 

A 7r (fc), a and 
The physical 



m 



zm - t X~ "t + Oim). Since A^fc) = 
7r have the same physical mass m 2 = m 2 . 
mass m 2 is obtained as a solution of 



'(27T) 



111 



2A 2 



(5.27) 



mt < 



In the region A ~ A cr , there is a nonzero solution 
2A 2 . 

In FIG. 4, we show the square masses of a and tt around 
A ~ A cr in the chiral and NJL limit (m — > 0, /i — > oo). In 
the broken phase (A > A cr ), m 2 = and m 2 = (2Af) 2 , 
while in the symmetric p hase (A < A cr ), m 2 = m 2 v is 
given as a solution to Eq. ( 5.27 ). The pion mass m 2 goes 
to zero in the limit A — * A cr — 0. 



we find the the effective pion-quark coupling = y 'Z„. 
Then, in the leading order of 1 /N expansion, f n is given 
by a fermion one-loop integral with appropriate Lorentz 
structure: 

ik »U = -{0\jg(0) \n(k)) 
d 4 p 



9*N / — ^ tr [r l5 S(p~ fc) 75 5(p)] 



(27T) 

2ig 7T NM k>*F E p l (k 2 



0,M, 



())• 



(5.24) 



her ml 




5 Fo r sm all but finite m, the condensate is evaluated from 
Eq. (|4.27|) as 



M = Mo + 



A C r 



1 



A2" 111 "SJ 



m + 0(m 



FIG. 4. Squared masses for m-+0 and /i — > oo scaled 
by A,^. 1 . Solid line is for cr and dashed line for 7r. For A > A cr 
we used the broken solution, whereas for A < A cr symmetric 
solution. 
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It is important to recognize that Eq. (5.27) implies the 
existence of tachyonic modes for A > A cr as we mentioned 
before. Indeed, if we assume A > A cr , we find a negative 
solution < 0. Therefore, if we choose the symmetric 
solution to the zero-mode constraint, then the resulting 
theory becomes unstable for A > A cr . So we must select 
the broken solution above the critical coupling. 



chiral charge Q^ F does not conserve. Since the r.h.s. is 
proportional to m 2 , it seems vanish in the chiral limit 
m 2 oc m — > 0. However, it does survive finite even in 
the chiral limit because the pion zero mode shows the 
singular behavior in the chiral limit: 



dx d 2 x±ir® oc — . 

m 



(5.31) 



C. Derivation of the PCAC relation and the 
nonconservation of Q\ ¥ 

As we discussed in Sec. III-B2, it is almost hopeless 
to treat the LF chiral current j£ f or t he massive case. 
Instead, we adopt the current j£ (3.7) which is much 
more tractable than j£ and gives the same null-plane 
charge. Then it is straightforward to derive the PCAC 
relation. Consider the divergence of jg in the u — > oo 
limit [see Eq. (|3.12fl ] 



« = 



„ N 
2m— 7r 

A 



2m- 



N 

T 



(5.28) 



where the normalized pion field 7r n was introduced so that 
its propagator be (fc) 



l/(k 2 



•mil 



If we use the pion decay constant (5.25) and the pion 



see Eq. (5.22)] 



mass (5.2C) [or equivalently m 2 = j^NZ w ], we obtain 
the PCAC relation, 

drfg = m*Uir u . (5.29) 
This is also consistent with Eq. Q5.24 ) [our normalization 

is <0|7T n |7r(fc)) = 1]. 

The important consequences of Sec. V-A were that i) 
the null-plane chiral charge Qg F is not conserved in the 
broken phase [<3g F , Hlf\ ^ 0, and that ii) the violation is 
proportional to the dynamical fermion mass M and does 
not disappear even in the chiral limit. Since we did not 
show explicitly the quantity because of its complexity, we 
here discuss it in more elegant way in the context of the 
PCAC relatio n. First of all, let us remember that the 
chiral charge ( |3.§| ) defined by jf? is equivalent to the one 
defined from j£ due to j§ — j§. This means that the 
time derivative of Q\ F should be the same even if we use 
the current j§ . Therefore, we ca n calc ulate the quantity 



d+Q\ F from the PCAC relation fl5.29| ): 



LF 



f^m 2 I dx d 2 x±ir n 



= f-nml J dx d' 



X±TT n 



(5.30) 



where 7r° is the zero mode of 7r n . If the current mass is 
not zero, the r.h.s. does not vanish in general and the 



Since f v i s pro portional to the fermion's dynamical mass 
Mo [Eq. ( 5.25 )1, we can confirm the result in Sec. V-A. 
That is, the null- plane chiral charge Q§ F does not con- 
serve d + Q\ F ^ even in the chiral limit. The singular 
behavior of Eq. ( 5.31 ) in the chiral limit has been pointed 



out by Tsujimaru and Yamawaki for spontaneously bro- 
ken theories p3| . They showed the necessity of intro- 
ducing the nonzero mass tting ^ for the NG boson in 
the broken phase and found the singular behavior of the 
NG-boson zero mode 6 cLS 



dx-d 2 x±4>o G 



NG 



We have confirmed their result for the dynamically bro- 
ken theory. 

Now let us verify the singular behavior ( 5.31 ). It is 



generally known that the zero-mode constraint for the 
NG boson becomes inconsistent in the broken phase un- 
less we introduce finite mass of the NG boson by han d 
as regularization |2S|| . Using the zero- mode mass ( 1.29 ), 
the zero-mode constraint for ttq is simply written as 



m 



ZM 



dx /i f(x , x±). 



(5.32) 



Suppose ttt-zm = and introduce the periodic 
boundary condition on ttq in the transverse direc- 
tions, then the transverse integral J d 2 x± of the 
zero-mode constraint leads to inconsistency^ 
j d 2 x± f^ L dx~ [i 2 f(x~ ,.xj_) 7^ 0, which suggests to in- 
troduce "zero-mode mass" mzM ^ 0. In our calculation, 
the origin of the the finite mass of the NG boson is the 
fermion's bare mass m [not the scalar mass fi/yX]. In- 
deed, m| M = n 2 mj\M in the broken phase. Therefore 
in the n — + oo limit, we find the singular behavior of the 
zero mode: 



1 



7TQ 



d 2 



dx n 2 f(x ,x±_) 



[ dx f(x ,arj_) 



(5.33) 



6 The mean-field result avoids this inconsistency due to 
J Mi'y^ m = 0. However, higher order calculation re- 
quires nonzero bare mass. 
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On the other hand, in the symmetric phase w| M = 
/i 2 (i — + 0{m) survives finite in the chiral limit. 
We called tuzm "zero-mode mass" , but it should not be 
confused with the physical pion mass m,. Both become 
nonzero due to nonzero bare mass m^O, but we have to 
calculate fermion's one loop to obtain the physical pion 
mass m-n- 

VI. CONCLUSION AND DISCUSSIONS 

We have studied a method of describing the dynamical 
chiral symmetry breaking on the LF. Our description is 
based on the idea in DLCQ that the symmetry break- 
ing is achieved by solving the "zero-mode constraints" , 
which already succeeded to some extent in describing the 
spontaneous symmetry breaking in simple scalar models. 
The point is that we can utilize this idea even for the 
dynamical symmetry breaking in fermionic systems if we 
introduce bosonic auxiliary fields for ^"5 and "J, and 
treat them as dynamical variables by adding their kinetic 
terms. Then the problem can be formulated such that we 
find a nontrivial solution to the "zero-mode constraint" . 
We exemplified this idea in the NJL model. The model 
we studied is the chiral Yukawa model, which reproduces 
the NJL model in the infinitely heavy mass limit of the 
scalars. Within this model, we showed in the massless 
case the equivalence between "light-front" chiral trans- 
formation and the usual one by classically solving the 
three (i.e. two zero-mode and one fermionic) constraints. 
This allowed us to construct the chiral current and 
charge Q5 F . Even if we solve the constraints classically, 
the resulting theory cannot have a symmetry breaking 
term. Quantum analysis showed that the zero-mode con- 
straint for a scalar boson a became the gap equation in 
nonperturbative treatment, which lead to nonzero con- 
densate {a) 7^ and equivalently to the chiral symmetry 
breakdown. We found the critical coupling A cr beyond 
which the fcrmion acquires nonzero dynamical mass. On 
the other hand, a perturbative solution could not give 
fermion condensate even in the quantum theory. 

The most important key of our description was the 
identification of the zero-mode constraint of a with the 
gap equation. This was suffered from a severe problem 
that the correct mass dependence disappears from the 
mode expansion. Of course this is a demerit of the LF 
formalism and we have to carefully incorporate mass de- 
pendence into e.g. ("J"]/) when we regularize its infrared 
divergence. It is suggestive that cutoff schemes with sym- 
metry consideration such as parity or rotational invari- 
ance lead to physically acceptable results. Contrary to 
such negative aspects, the mass-information loss has a 
useful and important aspect. It follows that the Fock vac- 
uum keeps invariant even if we change the value of mass. 
Therefore the vacuum does not change even though sym- 
metry breaking occurs and the fermion acquires dynam- 
ical mass. 



In our formalism, the vacuum is exactly the Fock vac- 
uum. The inclusion of dynamical scalar fields was nec- 
essary to clarify the structure of the Hilbert space and 
the triviality of the vacuum. The way of realizing bro- 
ken phase is that the vacuum is still trivial but the op- 
erator structure of the dependent variables changes. In 
other words, the "vacuum physics" in conventional for- 
mulation is converted into the hamiltonian through the 
dependent variables. The zero modes of scalars and the 
bad component ip_ are constrained variables and differ- 
ently expressed by physical variables depending on the 
phases. Related to this, the LF chiral transformation 
of the dependent variables also becomes unusual in the 
broken phase. Consequently, a seemingly contradiction 
between the triviality of the null-plane charge Q^ F \0) = 
and the chiral condensate (tyty) ^ is resolved. 

We further calculated masses of tt and a for both sym- 
metric and broken phases. In the broken phase, the mass 
of tt goes to zero in the chiral limit, which is consistent 
with the NG theorem. Our result is consistent with the 
GOR relation. If we substitute symmetric solution into 
the lagrangian, there appear tachyon modes for A > A cr . 
Therefore we can say that when A > A cr , physically real- 
ized phase is the broken phase. Certainly we have mass- 
less pion in the model, but it is very difficult to verify 
the NG theorem in general on the LF. This is because we 
have nonlocal interaction and because the chiral transfor- 
mation of the full fields explicitly depends on the model. 
Both of these arise from the fact that in the LF formal- 
ism the bad component of fermion and zero modes of 
scalars are constrained variables. However, we have suc- 
ceeded in deriving the PCAC relation. This was enabled 
by utilizing the current which was the Noether current 
defined for the massless fermion. The 'massive' current 
j'5 becomes complicated and it is almost hopeless to deal 
with it. The physics meaning of this discrepancy between 
jg and j'5 is still not clear but it seems that the usual 
current is favorable for discussing the "usual" chiral 
symmetry (not the "LF" chiral symmetry). 

One of the most important conclusion of our analysis 
is the nonconservation of the LF chiral charge Q5 F . This 
can be shown both by direct calculation of [QP^Hlf] 
using the unusual chiral transformation law and by uti- 
lizing the PCAC relation. In the broken phase, the chiral 
charge Q\ ¥ does not conserve d+Q$ F ^ even in the 
chiral limit. The singular behavior 1/m of the pion zero 
mode is essential to give a finite violation of d + Q\ F ^ 0. 

In our calculation, "nonperturbative" implied the 
mean-field approximation. This mean-field calculation is 
justified as the leading order approximation in the 1/N 
expansion. In principle, we can develop a systematic 1/N 
expansion to go beyond the mean-field result. Neverthe- 
less, the higher order will severely depend on the operator 
ordering and it is not clear whether the result with our 
specific ordering makes sense. If we want to go beyond 
the leading order, we have to determine the consistent 
operator ordering according to the criterion discussed in 
the text. Since this is a very difficult task in our model, 
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it should be examined in much simpler models such as 
1+1 dimensional Yukawa theory. 

It will be challenging to use other nonperturbative 
method to solve the zero-mode constraint. For example, 
Tamm-Dancoff approximation which truncates the Fock 
space into a few particle states will give some nontriv- 
ial results. Notice that our leading 1/N approximation 
corresponds to the two body truncation since multiquark 
states give higher order contribution. 

Our method here heavily relies on the introduction of 
auxiliary fields. So it seems natural to ask a question: 
Can we describe the chiral symmetry breaking without 
introducing the auxiliary fields? The answer is of course 
yes. Even though we do not have zero mode constraints, 
it is possible to describe the chiral symmetry breaking on 
the LF. Its explicit demonstration in the purely fermionic 
NJL model will be given in our next paper JL2J. As far 
as the NJL model is concerned, to solve the fermionic 
constraint becomes of great importance. 

There still remain many problems which cannot be dis- 
cussed in our model. One of them is the issue of renor- 
malization. In a renormalizable theory, if one introduces 
an infrared cutoff and excludes the zero mode degrees of 
freedom from the beginning, then the "vacuum physics" 
should be discussed as the problem of renormalization 
with nonperturbative infrared counter terms. Relation 
between such counter term approach 25 and the zero- 
mode approach presented here is not clear. We need 
further investigation for understanding how to describe 
the chiral symmetry breaking in LFQCD. 



V(a,w)/N= — (a 2 +7t 2 ) 
d 4 k 



+ i 



lndet$- + i757r)], (Al) 



which is the same result as that of the NJL model. There- 
fore evaluating the integral and differentiating V(a, n — 
0) with respect to c, we obtain the gap equation which 
determines the vacuum expectation value of a. 

Now let us turn to the mean-field approximation. Us- 
ing AB « (A)B + A(B) - (A) (B), the Yukawa interaction 
becomes 



Yukawa 
S3 (<t)#* + cr(**) 



(7r)*i7 5 * + 7r(*i7 5 *) + const., 



where (a) = 0(N°), a - (a) = 0{N- 1 / 2 ) and similarly 
for others. The leading order Euler-Lagrange equations 
for a and ir are (a) = -^(W) and (ir) = -^(^ij^). 
Therefore, in the mean-field approximation (^leading 
order of 1/N expansion), the chiral Yukawa model al- 
lows both fermion and scalar condensates. However, 
in higher order, fermion bilinears and scalars can in- 
dependently take their VEVs and the same relations 
do not necessarily hold. The Euler-Lagrange equation 
for fermion in the mean-field approximation becomes 
(i$ — m — (a) — (ir)i^)'i' = 0. Evaluating the fermion 
condensate in a self-consistent way, we obtain the gap 
equation. 
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APPENDIX A: THE CHIRAL YUKAWA MODEL 

Effective potential for scalars in the chiral Yukawa 
model (2.1) is easily calculated in the leading order of 
1/N expansion Exactly in the same way as in the 
Gross- Neveu model |2(J , we find that the leading contri- 
bution comes from the fermion one-loop diagrams. Note 
that the inclusion of kinetic terms for scalars has no ef- 
fect on the leading effective potential. Since the scalar 
propagator is 0(1/N), effects of the kinetic term (i.e. fi 
dependence) emerges from the next leading order. The 
effective potential V(a, ir) in the leading order is inde- 
pendent of /i and is given by 



APPENDIX B: CONVENTIONS 

We summarize our convention. We follow the Kogut- 
Soper convention First of all, the light-front coordi- 
nates are defined as 



1 



V2 



(x"±x i ), x\=x l (i = 1,2), (Bl) 



where we treat x + as "time" . The spatial coordinates x~ 
and xj_ are called the longitudinal and transverse direc- 
tions respectively. Derivatives in terms of x^ are defined 
by 



d+ = 



d 
dx±' 



For the 7 matrices, we also define 



1 



7 ± = -^(7°±7 3 ). 



(B2) 



(B3) 



It is useful to introduce projection operators A± defined 
by 



(B4) 



1G 



Indeed A± satisfy the projection properties A|. 



A±, A H 
projectors, 



A_ 



1, etc. Splitting the fermion by the 



(B5) 



we find that for any fermion on the LF, tp_ component 
is a dependent degree of freedom. tp + and tp- are called 
the "good component" and the "bad component" , respec- 
tively. 

In DLCQ, we set x~ finite x~ G [— L, L] with some 
boundary conditions on fields. Taking the periodic 
boundary condition, we can clearly separate a longitu- 
dinal zero mode from oscillating modes. The zero mode 
of some local function f(x) is defined by 



fo(x±) 



1 



dx f(x ,x±). 



The rest is the oscillating part: 

= f( x ) - fo(x±)- 

For some composite fields, we use the notation 
their zero modes: 

f(x)g(xj\ = t!t f dx-f(x)g(x). 
Jo 2LJ_ L 



(B6) 

(B7) 
]o for 

(B8) 



The inverse of the differential operator <9_ is defined as 

r-L 



1 f _1 

g~f( x ) = J L dy 2 €<yX ~ V ^ y ^ 
where e(x~) is a sign function 

r 1 (x- > 0) 
e(aT) = <^ (x- = 0) . 
{ -1 (x- < 0) 



(B9) 



(BIO) 



APPENDIX C: CONSTRAINTS AND THEIR 
CLASSICAL SOLUTIONS 

Since systems on the LF always have several con- 
straints, the LF quantization must be performed using 
Dirac's Hamiltonian formalism for constrained systems. 

The chiral Yukawa model has six primary constraints. 



Among them, consistency conditions for 9\ = n o 







and 82 — Una ~ (11$ is a conju gate momenta of gen- 



erate the zero-mode constraints ( |3.2| ), while 63 = H,p_ « 
generates the fermionic constraint (3.1). The consis- 
tency is calculated by 9\ = {9 i7 H^p + Xj9j}pB ~ with 
the canonical light-front Hamiltonian Plf = P~ 



P-= P F +P B +Py, 



(CI) 



P- = [ d 3 x {-* a i 7 -9_* Q - * a i 7 ± a ± * a + m* a * a } , 



(Fx 



N 



{(d^aod^oo + d ± ip a -d±<p a ) 



Pv 



d 3 x{(a + p ff )¥ tt tf° + (tt + Vir )tf rys* a } 



where J d?x = f_ L dx~ fd 2 x±. As a result, we find 
that no more constraints are generated from consistency 
conditions and that this system belongs to the second 
class. 

If we ignore ordering of the variables, it is not difficult 
to solve the constraints. First, the fermionic constraint 
jl) is solved as 



1 

id- 



(ij d± + m + a — 1^75) 7 + V^+ 



(C2) 



where dZ 1 is defined so that 4>-{ x ) also satisfies the an- 
tiperiodic boundary condition (see Appendix B). Note 
that however this solution still contains the zero modes 
00 and ttq and thus i s no t a complete solution. Substi- 
tuting ( |C2| ) into Eq. (|3.2|), we have equations only for er 
and 7Tq. Then the formal solution is given by 



&o( x ±)\ 
k (x ± ) J 

1 /i 2 



V2N 



V-\x x ) 



-(i"/ ± dj_ - m)V>+ 



,,at 1 



<y57r - «75<At 



^ + c.c. 



(C3) 



where the transverse differential operator D(xx) is 



T>(x±) 



EL 

A 



N ^2 



i9_ 



c f ) n 



The final expression for ^_ is reached after we insert 
Eq. ( |C3] ) into Eq. (C2). Though we have completely ig- 
nored the "ordering" in the classical treatment, the op- 
erator ordering becomes an issue in a quantum theory, 
which makes the analysis very complicated. 



APPENDIX D: PERTURB ATIVE SOLUTION TO 
THE ZERO-MODE CONSTRAINTS IN 
QUANTUM ANALYSIS 

Let us solve the zero-mode constraints in perturbation 
theory with the natural order ing in Eqs . (3.1) and ( p\2| ). 
In addition to the expansions ( 4.8 )-( 4.11 ) , it is convenient 
to define the expansion of "J: 



\ " 



A 



/ ] \Jjr«(*0 5 

t \ Ac 

n— 



(Dl) 
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where * a (°) = ip% + ipT> and * a ( fc ) = k) for k > 1. 

Knowing the lowest order solutions (4.13), we obtain 
the next order solution, 



r (1) 




$«(0) ( 1 ) ^a(O) 
»75 



/o(i) 1 1 Mi) ■ (i) A + / 

^ - 2 0L V° 75 J 7 ^ 



(D2) 
(D3) 



Similarly, if we know the solution up to n-th order, we 
easily obtain the (n + l)-th order solution because the 
constraint equation is written as follows 



2 / Jn+1) 
Jn+1) 



- f) 2 ( 



N 



a(Q) ( 1 



^a(O) 



V*75 



* Q(n) + h.. 



where ty a ( n ) for n ^ is 
1 1 



*° (n) ^^(-o n) -^ 5 )7 + ^ 



, (D4) 



(D5) 



and we have used i\r_ ^ip 1 ^ = for k, I ^ 0, etc. In 
this way, we can determine the solution order by order. 



APPENDIX E: NONLINEAR EQUATIONS FOR 
POLE MASSES 

With the Extended Parity Invariant (EPI) cutoff, the 
integral T(k 2 ,M 2 ) Eq. (5.18) is given as 



^ E pi(fc 2 ,M 2 ) 



(+) ^2K i x{\-x)-M A 

dx I dp\ 
Jo 



(-) 



7T 



1 



(2tt) 3 pl+M 2 - k 2 x{l - x) 
where integration limit xr±\ — (1 ± f3)/2 with j3 



(El) 



y/1 - 2M 2 /A 2 comes from 2A 2 x(l 



M 2 > 0. The 



integral is easily performed and the result is 
T E pi(k 2 ,M 2 ) 



1 (1 + 
8^ l0 HW 



4tt 2 



AM 2 - k 2 



k 2 



arctan (3 



k 2 



AM 2 - k 2 ' 



For k 2 ~ 0, we can approximate this as 
^ Pl( ^0,M*) = ^l„ s (I±|S 



4tt 2 



where (3q = ^/l — 2A/q / A 2 and we used 
linij^o x^ 1 arctan a; = 1. 

Using above, the nonlinear equations for the scalar 
masses are 



m 
AM 



m 2 ^ E pi(m 2 ,M 2 ) 
(ml -4M 2 )^ EP i(m 2 ,M 2 ) 



(E4) 



Note that when the chiral and heavy mass limit m — > 0, 
/i — * oo, we have solutions = and = 2M. 



(E2) 



/3b, (E3) [14 
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